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SUMMARY .  The  principal  equation  of  exterior  ballistics  has  a 
drag  term  which,  in  this  case,  is  proportional  to  the  square  of  the 
velocity  in  the  tangential  direction  of  the  projectile.  The  sensitivity 
coefficient  is  expressed  as  the  ratio  of  the  initial  elevation  angle 
deviation  to  the  initial  percentage  velocity  deviation.  The  work  in 
this  paper  is  to  find  analytically  the  sensitivity  coefficient  of  the 
exterior  ballistics  with  velocity  square  damping  which  comes  from  the 
nonlinear  air  resistance  for  a  projectile.  This  principal  equation  is 
integrated  analytically  in  obtaining  the  solution  for  tangential  veloc¬ 
ity  in  terms  of  the  elevation  angle,  together  with  all  the  necessary 
initial  conditions.  The  horizontal  range  and  the  vertical  range  are 
also  expressed  as  integrals  of  certain  function  of  the  elevation  angles. 
In  order  to  obtain  the  -sensitivity  coefficient  it  is  necessary  to  find 
the  perturbations  of  the  horizontal  and  vertical  ranges.  This  proce¬ 
dure  is  similar  to  that  of  evaluating  differentiation  under  the  integral 
sign.  The  perturbation  of  the  ranges  is  the  sum  of  the  perturbations 
due  to  the  initial  velocity,  the  initial  elevation  angle  and  the  impact 
elevation  angle.  By  setting  to  zeroes  the  range  perturbations  we  can 
group  the  coefficients  of  the  perturbations  into  two  separate  equations. 
The  ratio  of  the  perturbations  for  initial  elevation  angle  to  that  for 
initial  velocity  is  the  sensitivity  coefficient  for  exterior  ballistics 
that  we  are  seeking. 

I.  INTRODUCTION.  The  design  of  a  gun  involves  numerous  parameters. 
These  parameters  should  be  in  such  a  combination  that  the  best  first 
round  accuracy  is  given.  While  the  shell  leaves  the  gun  it  has  pertur¬ 
bations  for  the  muzzle  elevation  angle  and  the  muzzle  velocity.  The 
ratio  of  the  two  is  the  sensitivity  coefficient  of  the  interior  and  the 
exterior  ballistics.  It  is  desired  to  compensate  the  errors  due  to 
uncertain  changes  of  muzzle  velocity  by  the  automatic  response  of  the 
muzzle  elevation  angle  within  the  gun  system.  With  a  correct  design 
this  compensation  can  be  made  by  matching  the  exterior  ballistics  to 
the  interior  ballistics  through  the  analysis  of  gun  dynamics.  This 
process  is  called  passive  control  since  there  is  no  external  measure¬ 
ment  involved  nor  instrumentation  needed  for  control.  This  general 
problem  can  be  formulated  by  first  investigating  the  sensitivity 
coefficients  for  exterior  ballistics  with  velocity  square  damping. 
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II.  DYNAMICAL  EQUATIONS  FOR  TRAJECTORIES.  For  a  constant  mass 
travelling  in  a  vertical  plane  with  no  lift  and  applied  thrust,  but 
having  drag  and  velocity  vectors  contained  in  the  plane  of  symmetry  as 
shown  in  Figure  1,  the  dynamical  equations  of  motion  are  [1]: 

fl  Y 

■ii  -  vcos0  =  0  (1) 


-  vsin0  =  0 
dt 


m(gcos0  +  v  — )  =  0 
dt 

d2x  DcosO 


dt2  m 


where  m  = 

g  = 
D  = 
v  = 
0  = 
x  = 

y  = 


the  mass  of  the  projectile 

the  acceleration  due  to  gravity 

the  drag  of  the  projectile 

the  velocity  of  the  projectile 

the  path  inclination  (elevation  angle) 

the  horizontal  distance  of  the  projectile 

the  altitude  or  vertical  distance  of  the  projectile. 


(2) 

(3) 

(4) 


It  is  noticed  that  deviations  due  to  anomalies  in  the  azimuth  direction 
arenot  considered  here. 


By  differentiating  Equation  (1)  with  respect  to  t  one  obtains 


dzx 

dt2 


dt 


(VCOS0) 


Substituting  into  Eq.  (4)  we  have 

d  ,  DcosO 

—  (vcose)  =  -  -j- 

Solving  or  d0/dt  in  Eq.  (3)  one  obtains 

dQ  -gcosQ 
dt  v 


(5) 


(6) 


(7) 


Equation  (7)  indicates  that  the  differential  equations  can  be  transformed 
from  the  time  domain  in  t  to  the  angle  domain  in  0.  Equations  (6),  (l),and 
(2)  are  divided  by  Equation  (7)  in  achieving  this  transformation  as 


d(vcosQ)  _  £v 
d0  mg 


(8) 
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(9) 


dx  v2 
d0  "  "  g 


de 


—  tan0 
g 


(10) 


Equation  (8)  is  called  the  principal  equation  of  exterior  ballistics  [2]. 
It  can  be  integrated  if  the  drag  D  is  a  known  function  of  velocity  v. 


III.  TRANSFORMATION  OF  VARIABLES  UNDER  HEAD  WIND  DRAG.  The  head 
wind  drag  D  is  a  velocity  square  damping  term  given  as 


where 


D  =  mcv2 

c  =  cw(£  d2) (p/2) 


cw  =  the  dimensionless  resistant  coefficient 
d  =  the  diameter  of  the  projectile 
and  p  -  the  air  density. 


(11) 

(12) 


Thus  the  principal  equation  of  exterior  ballistics  (Equation  (8))  becomes 


_d_ 

d0 


(VCOS0) 


(13) 


A  further  transformation  of  the  dependent  variable  is  necessary  by  letting 

u  -  vcos0  (14) 


where  u  is  the  horizontal  component  of  the  projectile  velocity.  Then  the 
dynamical  Equations  (13),  (9)  and  (10)  become 


du 

d0 


£ 

g 


u3sec30 


dx 

d0 


sec20 


^  =  -  —  sec20tan0 . 
d0  g 


(15) 

(16) 
(17) 


To  simplify  further  the  form  of  the  dynamical  equations  another  trans¬ 
formation  of  the  independent  variable  is  made  by  letting 


q  =  tan0 


(18) 


where  q  is  the  projectile  slope. 
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Thus 


-3.  =  sec20  =  1+q2 
d0 

Equations  (15),  (16)  and  (17)  are  divided  by  Eq.  (19)  to  give 

du  =  —  u3  (l+q2f  ^dq 
g 

dx  =  -  —  dq 
:  g 

dy  =  -  —  qdq 
§ 


(19) 

(20) 
(21) 
(22) 


Solution  for  u  in  Equation  (20)  can  be  readily  integrated  in  closed  form. 
Solutions  for  x  and  y  can  be  expressed  in  the  form  of  integrals  once  u2 
is  obtainedo 

IV.  SOLUTION  FOR  HORIZONTAL  COMPONENT  OF  PROJECTILE  VELOCITY  AS 

FUNCTION  OF  TRAJECTORY  SLOPE.  The  solution  for  horizontal  component  of 

velocity  u  can  be  obtained  by  integrating  Equation  (20) 


-  j  (u‘2-V2)  =  {q(l+q2)1/2  +  «,n[q+(l+q2)1/2]} 


1/2.  -q 


2  g 


(23) 


where  qQ  equals  projectile  slope  initially  at  launch  and 


u  2  =  v  2sec-20  =  v  2(l+q  2)-1 
0  o  o  o  lo  J 


by  virtue  of  Equations  (14)  and  (19)  . 

Equation  (23)  can  be  written  as 

U-u02  |  [p (q)  -  P0(q0)} 


where 

and 


p (q)  =  q(i+q2)1/2  +  £n[q  +  (i+q2)1^2] 
PoCqo}  =  qo(1+qo2)1/2  +  *n[qo  +  C1+%2)1/2] 


Finally,  Equation  (25)  takes  the  form 


u2  = 


1+q, 


H(q,q0,vo%c/g) 

}2  +  1  -  H(q,qo,vo2,c/g)^ 


(24) 


(25) 

(26a) 

(26b) 

(27) 
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where 


Ep Cq)  -  P0(q0)] 


(28) 


H(q,q0,v02,c/g)  =  f  i~" 


It  is  noted  from  the  above  equation  that 


as 


For  the  case  with  no  air  resistance 


q  q„  H  -  0 
°  q=qo 


(29) 


as 


c  0  Hc=0  =  0  * 

which  implies  that  the  horizontal  component  of  projectile  velocity  u  at 
any  time  is  a  constant. 


(30) 


V.  SOLUTION  FOR  NONDIMENS IONAL  RANGE.  In  determining  the  range  x 
for  the  trajectory  the  closed  form  solution  of  u2  in  Equation  (27)  can 
be  substituted  into  Equation  (21)  to  obtain  the  solution  in  integral  fo 


Xi  '  X0  =  ' 


vo  r,  ,  fqi  H(q>q0>vo 

— — Jr  [(qi-q0^>  +  /  - dq 

g(l+qo  )  q  1  -  H 


(31) 


X  °  T  *  r 

x0  =  range  at  initial  point 
and  q^  =  projectile  slope  at  impact  ppint. 


where  x^  =  range  at  impact  point 
x~  =  ranee  at  initial  ooin 

.t  i: 

To  non-dimensionize  the  range.  Equation  (31)  is  divided  by  the  fact 
v02/g  as 

X^i'xo’vo)/A(qi’V>  =  Gx(WVQ2’c/g) 
where  the  nondimensional  range  is 

x(Xi,x0,v0)  =  (xi-x0)g/v02  , 

the  slope  function  is 

hr  -v  q0-qi 

A(q0>qp - r  » 

1+qo 

and  the  range  drag  function  due  to  air  resistance  is 

nr-  2  /  ^  ,  1  f^i  H(q.q0»V0  >C/§)  , 

S3  W  o  >c/^  *  1  -  Z  77  /  - T - 5 - dq  • 

q0-qx  q  i  -  h 

o 


or 


(32) 


(33) 


(34) 


(35) 
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It  is  noted  that  the  left  side  of  Eq.  (32)  contains  nc  drag  coefficient  c. 
It  is  ^heJtJ§.rip..,GXf.  that  i(s  a  function  of  H.,  which,  in.  turn  is  a  function  of 
the  drag  coefficient  "The  numerator  of  the  integral  in  Eq.  (35)  is  H. 
Since  Hc=q  =0  in  Eq.  (30),  the  range  drag  function  at  this  condition  is 


-  S  ■ 

rx? 


'  Gx(c=0)  =  i 


,,  separate  fopm  of  Jjq,.rfJ£2).  can  be  written  as 

gCxi-x0)  c1^2)  =  i  /u  _h_ 

v02  (.%-*&  "  q  1-H  q 


(36) 

5  >  -_7 

(37) 


■'  VI 0  -VARIATION  OF  THE  NOND I  MEN  S I ONAL  RANGE  AND  THE  SLOPE  FUNCTION. 

In  order  to  obtain  a  first  round"  hit  of  the  target  one  of  the  conditions 
is  that  the  variation  of  the  range  should  be  zero,  i.e.,  from  Eq.  (33) 

rf,f  ?  •.  ’>•  _  , :  .  •»  r  ’  ,  ,j  .  4  - 

6(x^-x0)  =  0  .  (38) 

We  take  the  perturbation  for  the  nondimensional  range  from  Eq.  (33) 

•  " 6X  _  ($Cxi-xo)  26v0 

X 

.in  l  r:  -Wi  ■•!  -  ■ 


as 


xi"xo 


6X  _  „  26vo  .«.  ' 

X  “  ”  v  * 

»-  (  I  y  a  '»*<;. ;  i  r-p  ;;;  •  >  <i(|M  t>  ('  *  *  ;  . 

The  variation  of  the  slope  function  A  in  Eq.  (34)  becomes 


6A 

X 


6qi  6q0  2q06q0 

-  + 


q  -q. 

no  ni 


q  -q .  1+q 

no  ni  c 


(39) 


(40) 


Next,  taking  the  variation  of  Eq.  (32)  and  using  the  expressions  given 
in  Eq.  (39)  and  (40)  we  have  '  '  •  > 

6Gy 


or 


r 


6X 

6A 

•  X 

A  ' 

26v0  ,  6qi 

:  Ho 

\ ■  \  v  f  X.— 

y^o:  .  qo"qi 

qo"qi 

+  ^ 


1+qo  V 


(41) 

(42) 
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It  is  noticed  that  in  the  absence  of  air  damping  the  range  drag 
function  is  unity  from  Eq.  (36)  and  the  variation  <5GX  =  0.  Uqder  this 
conditiqn  the  solutions  for  Eq.  (42)  were  given  by  the  author  in  the  paper 
entitled,  "On  the  Sensitivity  Coefficient  of  Exterior  Ballistics  and  Its 
Potential  Matching  to  Interior  Ballistics  Sensitivity".  This  paper  was 
presented  at  the  Second  U.S.  Army  Symposium  on  Gun  Dynamics,  September 
1978.  With  the  velocity  square  damping  the  variation  of  Gx  is  not  zero, 
i.e. , 


«GX  f  0. 


(43) 


YjL_  VARIATION  OF  THE  RANGE  DRAG  INTEGRAND.  The  range  drag  fqnction 
in  Eq. (35)  can  be  written  as 

l  ^i 

Gx(qi’Vvo2’c/S)  =  1_  T-'-T-  f  vV,c/g)dq  (44) 

4°"qi  % 

where  the  range  integrand  is 


F(q*q0,v02,c/g) 


H(q,q0,Vo2,c/g) 
1  -  H 


(45) 


The  variation  of  Gx  involves  the  initial  velocity  v0,  the  initial  slope 
q0,  and  the  impact  slope  qj.  r 

Eq.  (44)  has  the  parameters  q0  and  q^  in  the  denominator  as  well  as 
in  the  integral  of  F,  By  chain  rule  we  have 


6Gx  =  "  a- -a"  1  F(q»q0>vo2>c/g)dq] 

4°  4i  % 

,  fa 

-(-i)(q0-qi)_  6(q0-qi)/  F(qJq0»v02>c/g)dq  .  (46) 

qo 

The  variation  of  the  integral  of  F  is  given  in  the  next  section. 

VIII.  VARIATION  OF  THE  RANGE  DRAG  INTEGRAL.  The  parameters  in  the 
integral  are  q0,  qj  and  v0.  The  variation  of  the  range  drag  integral 
follows  the  rules  of  differentiation  under  the  integral  sign. 


/ 


8 


q . 

fi[J  1  F(q»q0»vo2»c/s)dq] 
q„ 


=  [/  1  dq]6v0  +  F(q=qi,q0,v02,c/g)6qi 


q  8vo 
n0 


[/. 


i  9F_ 

9q0 


dq]6q0  -  F(q=qo,qo>V02,c/g)6qo 


The  last  term  of  this  equation  is  zero  by  virtue  of  Equation: 
(45).  Substituting  Eq.  (47)  into  Eq.  (46)  gives 


SG  -  [- 


1 


q0-qi 


qi  ^L 

9vr 


dq]  5v 


+  [ - —  F(q=q  )  -  - - l—j  J1  F  dq]6q. 

.  ■  q0-qi  (%-qi)  q 


+  [- 


i 
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li  3F 


q0-qi  q  ^o 
n0 


dq  + 


To 

A 


(q0-qi)  q 


/  F  dq] 6q 


IX.  EVALUATION  OF  THE  PARTIAL  DERIVATIVES  OF  THE  RANGE 

The  partial  derivatives  of  F  with  respect  to  v0  can  be  found 
(45)  and  (28). 

9H  9H  9H 


BF  (1-H)  9v0  -  H(-  9v0) 
9v„ 


(1-H)' 


<jv£ 

(1-H)2 


where 


9H0  c  2v0 


1  3v0  g  l+q0 
Combining  the  above  we  have 

9F 


t  [p(q)-p(q0)l  =  ~  h 


H 


9v0  v0  (1-H) : 


(47) 

(29)  and 


(48) 

DRAG  INTEGRAND., 
by  using  Eqs. 

(49) 

(50a) 

(50b) 
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Similarly  the  partial  derivatives  of  F  with  respect  to  qQ  is 


3F 

3q0 


where 


9H_ 

9% 


9H 

(1"H) ' 


■  7  VC^hi^o2)  2%[p(q)-p0cq0)] 


c 

+  — 


v-  dP0 


8  i-*q02 
2q0 

H  - 


(-1) 


dq0 

»2  dp0 


iV  g  lHq02  dq° 

Substituting  Eq.  (52)  into  Eqe  (51)  one  obtains 


3F  2q0 


9q0  i+q02  (i-H)2 


dPo 


g  (l+q02)  dq0  (1-H) 


(51) 


(52) 


(53) 


X.  VARIATION  OF  THE -RANGE  DRAG  FUNCTION.  By  substituting  Eqs.  (45), 
(50)  and  (53)  into  Eq.  (48),  we  have 


6Gx 


=  [- 


1  M  H 

qG-qi  q0  (1-H)2 


dq] 


26v0 


2!sai - ■ —  rqi  jl  dq)  Jsi_ 

i_^q=q^  q0-qi  q0  qo_qi 


+ 


[■ 


2q0  jqi 


H 


i+q0' 


(1-H) 


dq 


/ 


qi 


H 


q0-qi  q 


dq 


+  £  vo2  dPo  .qi  dq  j  <Sqo 
g  1-flo  2  dqo  q  (!-H)2  q0-qi 


(54) 


It  is  noted  that  the  difference  of  the  end  slope  is  not  zero,  i.e., 
qQ-q^  ^  0.  Therefore,  the  problem  does  not  become  singular.  We  have 
expressed  the  variation  6GX  in  terms  of  the  variational  parameter  6v0, 
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and  6q0p  However,  the  variational  parameter  6q^  at  the  impact  point 
is  not  known  explicitly  and  must  be  eliminated  by  using  another  variation 
in  the  direction  of  the  elevation  y. 


XI.  VARIATIONAL  EQUATION  F,0R  THE  RANGE.  By  substituting  Eq.  (54) 
into  Eq.  (42)  and  grouping  the  coefficients  for  the  variational  terms, 
we  have 


where 


and 


i  .  i 


<Sqo 

v0  qi  q0-qi  *q°  q0-qi 


iq.  =  -i  - 


1 

X  Hi 

1  1  T  ^  1  Hi 

Gx  %-ii  ll(,0’qi)  '  cx  wq- 


=  1  *  ^  ~V .’ufVii)  *  t 


Gx  ^o’^i 


1  ,  2q0 

*  Hq0" 


(55) 

(56) 

(57) 


1  c  vo2  dPo  2q0 

K  WVV  -  1^7  tq^qp 


(58) 


In  turn,  the  integral  1^,  I-^ 
as  follows. 

q. 

Ill(qo,qi)  =  / 

^O 

112^d^0 

I02(qo‘qP 

a  Hi  = 

and 

Gx  -  1 


, and  Iq2»  and  other  terms  are  given 


H(q,q0,v0  ,c/g)  ^ 

(59) 

1-H 

=  jqi  H  T  a,  , 
q  U-H)2 

Mo 

(60) 

;qi  1  a 

~  I  ,  dq  » 

(61) 

q0  U-H)2 

H(q=qi) 

(62) 

- —  T11  • 

qo-qi  il 

(63) 
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It  is  noted  that  6q^  in  Eq.  (55)  has  to  be  eliminated  in  solving  the 
sensitivity  problem.  Similar  variation  equation  may  be  obtained  by 
considering  the  variation  of  the  elevation. 


XII.  THE  SOLUTION  FOR  ELEVATION.  The  differential  equation  for 
elevation  was  given  in  Eq.  (22)  and  the  solution  for  u  is  in  Eq.  (27). 

Substituting  Eq.  (27)  into  Eq.  (22)  gives 


dy 


g(i+q„  ) 


ri  H(q,q0,v0  ,c/g) 

[1  +  - - - - - -  ]  qdq 

1  -  H 


Integrating  the  above  one  obtains 


g(i+q0  ) 


l-H 


(64) 


(65) 


Rearranging  yields  the  relationship  between  the  range  Y,  the  end  slope 
function  A,  and  the  elevation  drag  function  Gy. 


Y(yi,y0,vQ)/A(qi,q0)  =  J  (q0+qi)  +  Gy(qi,q0,v02,c/g) 
where  the  nondimensional  elevation  is 

g(yi-yo) 


Y(yi,y0>vo)  = 


vr 


(66) 


(67) 


A  is  given  in  Eq.  (20)  and  the  elevation  drag  function  is 


1  qi 

Gy(Wvo2,c/g)  =  ■  /  qp Cq,q0 ,v02 ,c/g)dq 

0  1  % 

A  separate  form  of  Eq.  (66)  may  be  written  as 

g(yi-y0)  #  q0-qi  1  1  ,q^ 

t— C  =  7  (q°*qi)  *  ^7  Jq  qF  dq 

O 


(68) 


(69) 


Note  that  the  left  side  of  Eq„  (66)  contains  no  drag  coefficient  c. 

It  is  the  term  Gy  which  is  a  function  of  drag  coefficient  c. 

For  c  =  0  Gy  -  0  •  (70) 
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XIII.  TERRAIN  SLOPE  FROM  LAUNCH  POINT  TO  TARGET  POINT.  If  Eq.  (69) 
is  divided  by  Eq„  (37)  with  the  aid  of  Eqs.  (35)  and  (68),  one  obtains 

Yi-Xo  .  !(l/2Xqo+qi)+Gy 

- A  m  =  - -  (71) 

xi“xo  =  Gx 

where  m  is  the  terrain  slope  from  launch  point  to  target  point,  a  constant 
parameter o  Therefore, 


(l/2Xq0+qi)  +  Gy  =  mGx  .  (72) 

It  is  noted  that  for  m  =  0, 

qi  =  -q0  -  2Gy  •  C73) 

‘From  Equations  (36),  (70)  and  (72)  we  have  for  c  =  0, 

^i  +  q0  ■  m  •  (74) 

••  '•  (  V  1 

We  use  Equation  (71)  to  find  the  variational  equation  for  the  elevation. 
Taking  the  variation  of  Eq.  (72)  for  any  given  m,  we  have 

(l/2X«q0+6qi)  +  <5G  -  mG  -  0  (75) 

Gx 

where  6GX/GX  is  given  in  Equation  (42)  and 

6Gy  =  “  trir  6[/qi  qF  dq] 

4°  4i  % 

qi 

-  (-^(qo-qp'^Cqo-qi)/  qF  dq  (76) 

qo 

is  obtained  from  Eq.  (68). 

It  is  noted  that  for  c  =  0,  both  6Gy  and  6GX  are  zero  in  Eq.  (75). 
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It  can  also  be  proved  that  the  result  for  6Gy  is 


<5G, 


[- 


1  fqi  qH 
qQ-qi  q  (1"H) 1 


dq] 


[- 


qiHq=qi 

1-Hq=q. 


—  ^  dq] 

-q,  1-H  HJ 


q0-qi 


2  6vr 


foii 

q0-qi 


♦  ,Z3£_  j"1  _SH 


q0 


j  1  r  1  H  j 

■5-  dq  +  -  I  -  dq 

(1-H)2  qQ-qi  J  1-H 


+  £  vo2  ^Pq  qdq  -»  ^0 

g  ^q02  dqQ  q  (1-H) 2  qQ-qi 

*  O 


(77) 


XIV.  VARIATIONAL  EQUATION  FOR  THE  ELEVATION.  By  substituting  Equa¬ 
tions  (77)  and  (42)  into  Equation  (75)  and  grouping  the  coefficients  for 
the  variational  terms,  we  have 


where 


and 


26v0 

J„  -  + 


6qi 


V  vr 


Jqi  qo-qi  +  Jq°  q0-qi 


=  0 


Jv  “  ” 


qn-q 


Jl2  +  mGx 


O  nl 


Jqi  -  -  \  cq0-qi)  -  -1- J 


qiHi 


q0-qi  11  1-% 

2q0 


-  mG^ 


Jq° =  7  (qo-qi^ +  ^7  Jn  +  J 


i+V  12 


(78) 

(79) 

(80) 


C  v02  dPo 
+ - 5- - J 

s  i+q0  dq0 


02 


_  r  ,  2q0(q0-qj) 

mG  [-1  +  - - - 5 — 

x  i+q0 


(81) 


In  turn,  the  integral  J^,  J^2  and  d02»  and  °ther  terms  are  given  as 
follows. 

q. 

JH  =  /  1  qH(q,q0,VQ2,c/g)/(l-H)  dq  , 


(82) 
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(83) 


and 


(84) 


i 


G 


x 


(85) 


XV.  THE  SENSITIVITY  COEFFICIENT  FOR  EXTERIOR  BALLISTICS.  Eliminating 
6qi/(qj^q^)£rom  Equations  (55)  and  (78)  we  have 


From  which  one  obtains  the  sensitivity  coefficient  through  the  aid  of 
Equation  (19) 


39q  _  (ly/lqj)  ~  (Jy/Jqj)  ^-2(q0-q^)^ 
6v0/vo  ’  (Iq0/Iqi)  ~  (Jq0/Jcii)  1+q02 


s 


(87) 


It  is  noted  that  Equation  (87)  requires  the  evaluation  of  the  integrals 
I  and  J,  which  are  given  in  Equations  (56)  through  (63)  and  Equations 
(78)  through  (85). 

XVI,  CONCLUSIONS.  The  following  results  are  concluded  in  this  paper: 

1.  The  principal  equation  of  exterior  ballistics  is  derived  with 
the  trajectory  slope  as  independent  variables. 

2.  The  closed  form  solution  for  the  horizontal  component  of 
trajectory  velocity  is  determined  for  the  case  of  exterior  ballistics 
with  velocity  square  damping. 

3.  The  nondimensional  range  is  obtained  in  terms  of  an  end  slope 
function  and  a  range  drag  function. 

4.  Variations  of  the  nondimensional  range  are  expressed  as  variations 
of  launch  velocity. 

5.  Variations  of  the  range  drag  function  are  in  terms  of  the  varia¬ 
tions  of  the  range  drag  integral. 
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6.  The  range  drag  integral  has  parameters  in  the  integrand  as 
well  as  the  upper  and  lower  limits.  The  variations  of  this  integral 
are  found. 

7.  The  partial  derivatives  of  the  range  drag  integrand  are  eval¬ 
uated. 

8.  The  variational  equation  for  the  range  is  in  terms  of  elements 
involving  three  integrals  as  coefficients  of  three  variational  para¬ 
meters. 

9.  The  variational  parameters  are  that  of  launch  velocity,  the 
launch  elevation  angle,  and  the  impact  elevation  angle. 

10.  The  average  of  the  end  slopes  is  equal  to  the  terrain  slope 
times  the  range  drag  function  minus  the  elevation  drag  function. 

11.  Variations  of  the  nondimens ional  elevation  are  expressed  as 
variations  of  the  end  slopes  and  the  variations  of  the  drag  function. 

12.  The  variational  equations  for  the  elevation  are  determined 
similar  to  that  for  the  range. 

13.  Eliminating  the  variations  of  impact  slope,  fiq^,  from  the  set 
of  two  variational  equations  gives  the  ratio  of  the  coefficients  of 
6v0/v0  and  6q0/(q0-qi)» 

14.  The  sensitivity  60o/(6vo/vo)  may  be  obtained  by  dividing  this 
ratio  6q0/ (6v0/v0)  by  the  quantity  (l+qD2). 

Numerical  calculations  of  this  problem  will  be  carried  out  in  the 
future . 
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